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Abstract: We study Nf D6-brane probes in the supergravity background dual to D4- 
branes compactified on a circle with supersymmetry-breaking boundary conditions. In the 
limit in which the resulting Kaluza-Klein modes decouple, the gauge theory reduces to non- 
super symmetric, four-dimensional QCD with colours and A^f <C A^c flavours. As expected, 
this decoupling is not fully realised within the supergravity/Born-Infeld approximation. For 
A'^f = 1 and massless quarks, = 0, we exhibit spontaneous chiral symmetry breaking by 
a quark condensate, (V'V') 7^ 0) fi^id the associated massless 'pion' in the spectrum. 
The latter becomes massive for > 0, obeying the Gell-Mann-Oakes-Renner relation: 
= —m^{'ijjip) / f^. In the case A'^f > 1 we provide a holographic version of the Vafa- 
Witten theorem, which states that the U{Ni) flavour symmetry cannot be spontaneously 
broken. Further, we find — 1 unexpectedly light pseudo-scalar mesons in the spectrum. 
We argue that these are not (pseudo-)Goldstone bosons and speculate on the string mechanism 
responsible for their lightness. We then study the theory at finite temperature and exhibit 
a phase transition associated with a discontinuity in (ipip) (T) . D6/D6 pairs are also briefly 
discussed. 
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1. Introduction and summary of results 

The AdS5/CFT4 correspondence relates string theory on the near- horizon background of 
D3-branes to the conformal field theory living on their worldvolume [1] . Inspired by this corre- 
spondence, Witten proposed a construction of the holographic dual of fom-dimcnsional, pure 
SU(N^) Yang-Mills (YM) theory [2]. One starts with A''^. D4-branes compactified on a circle 
of radius M^k, and further imposes anti-periodic boundary conditions for the worldvolume 
fermions on this circle. Before compactification, the D4-brane theory is a five-dimensional, 
supersymmetric SU{N^) gauge theory whose field content includes fermions and scalars in 
the adjoint representation of SU{N^), in addition to the gauge fields. At energies much 
below the compactification scale, Mkk, the theory is effectively four-dimensional. The anti- 
periodic boundary conditions break all of the supersymmetries and give a tree-level mass 
to the fermions, while the scalars also acquire a mass through one loop-effects. Thus at 
sufficiently low energies the dynamics is that of four-dimensional, massless gluons. 

The D4-brane system above has a dual description in terms of string theory in the near- 
horizon region of the associated (non-supersymmetric) supergravity background. Unfortu- 
nately, as observed in [2], the Kaluza-Klein (KK) modes on the D4-brane do not decouple 
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within the supergravity approximation. For example, the mass of the hghtest gluebaU is 
of the same order as the strong couphng scale. In this sense, there is no energy region, 
Aqcd ^ £' ^ -^KK, in which a description in terms of weakly-coupled gluons is appropriate. 
Nevertheless, the qualitative features of the glueball spectrum agree with lattice calculations 
[3], and the supergravity description suffices, for example, to exhibit an area law for Wilson 
loops [2]. 

Since most of our phenomenological understanding of low-energy QCD comes from the 
study of mesons and baryons, it is interesting to construct a holographic dual of a gauge 
theory whose low-energy degrees of freedom consist not only of gluons but also of fundamental 
quarks. In the context of AdS/CFT, fundamental matter can be added to the gauge theory 
by introducing D-brane probes in the dual supergravity background [4, 5, 6, 7, 8, 9, 10, 11]. 
This has been recently exploited to study mesons holographically in several examples of 
gauge/gravity duals [12, 13, 14, 15]. The goal of this paper is to study a gauge/gravity dual 
in which the gauge theory reduces to non-supersymmetric, four-dimensional QCD in the limit 
in which the KK modes decouple. The construction is as follows. 

Consider the D4/D6 system with the branes oriented as described by the following array: 

A^e D4: 1 2 3 4 

TVf D6: 1 2 3 . 5 6 7 

Note that the D4- and the D6-branes may be separated from each other along the 89- 
directions. This system is T-dual to the D3/D5 intersection, and in the decoupling limit 
for the D4-branes [16] provides a non-conformal version of the AdS/dCFT correspondence 
[5, 17]. On the gauge theory side one has a super symmetric, five-dimensional SU{N^) gauge 
theory coupled to a four-dimensional defect. The entire system is invariant under eight super- 
charges, that is. A/" = 2 supersymmetry in four-dimensional language. The degrees of freedom 
localized on the defect are Nf hypermultiplets in the fundamental representation of SU{N^), 
which arise from the open strings connecting the D4- and the D6-branes. Each hypermultiplet 
consists of two Weyl fermions of opposite chiralities, V'l and tp^,, and two complex scalars. 

As discussed above, identifying the 4-direction with period and with anti- 

periodic boundary conditions for the D4-brane fermions, breaks all of the super symmetries 
and renders the theory effectively four-dimensional at energies E ^ Mkk- Further, the ad- 
joint fermions and scalars become massive. Now, the bare mass of each hypermultiplet, m^, 
is proportional to the distance between the corresponding D6-brane and the D4-brancs. Even 
if these bare masses arc zero, we expect loop efi'ects to induce a mass for the scalars in the 
fundamental representation. Generation of a mass for the fundamental fermions is, how- 
ever, forbidden by a chiral, U{1)a symmetry that rotates V'l and ipn with opposite phases. ""^ 
Therefore, at low energies, we expect to be left with a four-dimensional SU{N^) gauge theory 
coupled to Nf flavours of fundamental quark. 

In the dual string theory description, the D4-branes are again replaced by their supergrav- 
ity background. In the so-called 'probe limit', Nf ^ N^, the backrcaction of the D6-brancs on 
^This symmetry is broken by instanton effects, but these vanish in the large-A'c limit in which we work. 
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this background is negligible and hence they can be treated as probes. The D6-brane world- 
volume fields are dual to gauge-invariant field theory operators constructed with at least two 
hypermultiplet fields, that is, meson-like operators; of particular importance here will be the 
quark bilinear operator, tpijj = %bitp^, where ■^^ = + V'r ^^^^ ^ = 1, . . . , A^f is the fiavour 
index. In the string description the U{1)a symmetry is nothing but the rotation symmetry 
in the 89-plane. In addition to acting on the fermions as explained above, this symmetry also 
acts on the adjoint scalar X = + iX^ by a phase rotation. 

Having presented the general construction, we now summarise our main results. 

We begin in section 2 by considering the case of a single D6-brane. For = 0, we expect 
the chiral U{1)a symmetry to be spontaneously broken by a chiral condensate, {'4'^p) ^ 0, and 
there to be an associated massless, pseudo-scalar Goldstone boson in the spectrum, as in 
QCD with one massless flavour. In QCD, this is the r]' (which becomes massless in the 
large- A^c limit), but, in an abuse of language, we will refer to it as a 'pion'. We are indeed 
able to show that the string description provides a holographic realisation of this physics, 
by showing that the U{1)a symmetry is spontaneously broken by the brane embedding, as 
in [14]. Moreover, we numerically determine the chiral condensate for an arbitrary quark 
mass. We find {■ipip){m^ = 0) / 0, as expected, and {tpip) (m^) oc l/rWq for — > oo, again in 
agreement with field theory expectations (which we briefly review). 

In sections 3 and 4 we study the meson spectrum. We first show analytically, in section 
3, that at mq = there is exactly one normalisable, massless pseudo-scalar (in the four- 
dimcnsional sense) fluctuation of the D6-brane. This open-string mode corresponds precisely 
to the Goldstone mode of the D6-brane embedding. If mq > 0, this pion becomes a pseudo- 
Goldstone boson. We are able to show, analytically, that its squared mass scales linearly 
with the quark mass, oc m^, in the limits of small (section 4) and large (section 3) 
quark mass. With numerical calculations we then confirm that, in fact, such a linear relation 
holds for all quark masses. For small mq this linear relation is in perfect agreement with the 
Gell-Mann-Oakes-Renner (GMOR) relation [18], 

^,^_m^^ (1.2) 

which gives the first term in the expansion of the pion mass around mq = 0. In section 4, we 
compute the chiral condensate and the pion decay constant analytically (at mq = 0) , and with 
these results we arc able to verify that, for small quark mass, the mass of oTir pion precisely 
satisfies the GMOR relation. In the opposite limit, mq — oo, we show in section 3 that 

^7r = «-F^' (1-3) 

where a is a pure number; for the lightest mesons, we find a ~ 20. 

The remaining mesons (D6-brane excitations) are studied in section 3. For mq < M^^^y, 
where 

M.u.y = gl^KM^^ , (1.4) 
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they all have masses of order the compactification scale, Mkk, and so the (pseudo-)Goldstone 
boson dominates the infrared physics in the regime of small mg. Of course, this result signals 
the lack of decoupling between the KK and the QCD scales within the supergravity/Born- 
Infeld approximation, as expected. For > -^susyj supersymmctry is approximately re- 
stored, the spectrum exhibits the corresponding degeneracy, and all meson masses obey a 
formula like (1.3) with appropriate values of a. The fact that supersymmctry is restored at 
mq ~ Mgusy suggests that the effective mass of some of the microscopic degrees of freedom 
must be at least of order M^^^y. We close section 3 with an examination of certain stability 
issues for the D6-brane embeddings. 

In section 5, for the case of multiple flavours, A^f > 1, we provide a holographic version 
of the Vafa-Witten theorem [19], which states that the U{Nf) flavour symmetry cannot be 
spontaneously broken if > 0. In the holographic description this is realised by the fact 
that the Nf D6-branes must be coincident in order to minimise their energy. 

In section 6, we examine the theory at finite temperature. If > M^^^y, the theory 
exhibits two phase transitions: the first one is the well-known confinement /deconfinement 
transition at T = T^ecoji! ~ -^kk [2] , and the second one is a phase transition at 



characterised by a discontinuity in the chiral condensate, {ipilj){T), and in the specific heat. 
If TTiq < Mgusy the second transition does not occur as a separate phase transition. 

In section 7 we discuss six-brane embeddings that, from the brane-construction viewpoint, 
correspond not to D4/D6 intersections but to intersections of D4-branes with D6/D6 pairs. 
From the gauge theory viewpoint these configurations correspond to having a defect/anti- 
defect pair. 

We close with a brief discussion of our results in section 8. In particular, one result of the 
analysis of the TVf > 1 case in section 5 is the existence of massless, pseudo-scalar mesons 
in the gauge theory spectrum if rriq = 0. (We expect the inclusion of sub-leading effects in 
the expansion to generate small masses for these particles.) Although it is tempting 

to interpret these particles as Goldstone bosons of a putative, spontaneously-broken U {Nf)p^ 
chiral symmetry, we argue that this interpretation cannot be correct: even at large N^, only 
one of them is a true Goldstone boson. Thus we are unable to find any obvious reasons, 
from the viewpoint of the four- or five-dimensional gauge theory, for these scalars to be light. 
We speculate that, from the string viewpoint, ten-dimensional U{N;) gauge invariance (as 
opposed to just seven-dimensional gauge invariance on the D6-branes) is responsible for their 
lightness. 

2. Chiral symmetry breaking from D6-brane embeddings 

In this section we will show how the spontaneous breaking of the U{1)a chiral symmetry 
expected from the field theory side is realized in the string description. 




(1.5) 
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2.1 The D4-soliton background 

The type IIA supergravity background dual to D4-branes compactified on a circle with 
anti-periodic boundary conditions for the fermions takes the form 

ds' = [^j {v,. dx^dx'^ + f{U)dr^) + (^^ j + R'/^U^^ d^l , (2.1) 

^'' = 9s(^] , F, = '^e,, fiU) = l-^. (2.2) 



Rj ' ^ V4 ^ ' ' C/3 • 

The coordinates = {x^, . . . ,x^} parametrize the four non-compact directions along the 
D4-branes, as in (1.1), whereas r parametrizes the circular 4-direction on which the branes 
are compactified. dil^ ^nd are the S'0(5)-invariant line element and volume form on a 
unit four-sphere, respectively, and = Stt^/S is its volume. U has dimensions of length and 
may be thought of as a radial coordinate in the 56789-directions transverse to the D4-branes. 
To avoid a conical singularity at ?7 = C/kk, t must be identified with period 

47r i?3/2 

Following the nomenclature of [21], we refer to this solution as the D4-soliton. 

This supergravity solution above is regular everywhere and is completely specified by 
the string coupling constant, Qs, the Ramond-Ramond flux quantum (i.e., the number of 
D4-branes), A'^,, and the constant J7kk- The remaining parameter, R, is given in terms of 
these quantities and the string length, £5, by 

R^ = TrgsNJ^,. (2.4) 

If UjiK is set to zero, the solution (2.1, 2.2) reduces to the extremal, 1/2-supersymmetric 
D4-brane solution. Hence we may say that [/kk characterises the deviation of the D4-soliton 
from extremality. 

The SU{N^) field theory dual to (2.1, 2.2) is defined by the compactification scale, Mkk, 
below which the theory is effectively four-dimensional, and the four-dimensional coupling 
constant at the compactification scale, gym- These are related to the string parameters by 

o 7-rl/2 o 7-/1/2 /or/ \ V2 

KK = ttt;^ = 77^: .1/0 ' fYM = 3Vvr -^-^ . (2.5) 



2 i23/2 20F (5,iVjV2^3/2 ' .YM V ^ 

The first equation follows from the fact that r is directly identified with the compact direction 
in the gauge theory, so Mkk = 27r/(5T. The second equation follows from the fact that 
= (27r)^(?s^s and the relation between the four- and five-dimensional coupling constants, 
gr^^j = g'^^/Sr, which can be seen by expanding the Born-Infeld action of the D4-brane and 
adopting the standard normalization for the gauge field kinetic term, —F'^/Ag'^. 
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It is useful to invert tliese relations to express the string parameters in terms of the gauge 
theory ones: 

^ - ' - 2^M^s ' " 9^--^^ ^^^^^ ■ ^^-^^ 

The string length will cancel in any calculation of a physical quantity in the field theory. For 
example, the QCD string tension is^ 

Now we may ask in what situation the D4-soliton provides a reliable background in which one 
can study the dual field theory using only classical supergravity. First, we must require that 
the curvature is everywhere small compared to the fundamental string tension. This ensures 
that higher derivative string corrections to the low energy equations of motion are negligible. 
The maximum curvature in the geometry given by (2.1) occurs at precisely U = C/kk; where 
the curvatures are of the order (?7kk-R^^)~^''^- Hence we require 

rrV2 o3/2 

^^^^ gl^N^ » 1 , (2.8) 



s 



where we have applied the results in eq. (2.6). Therefore the restriction to small curvatures 
corresponds to a large 't Hooft coupling in the effective four-dimensional gauge theory, pre- 
cisely as in the conventional AdS/CFT correspondence [1]. Further, to suppress string loop 
effects, we must also require that the local string coupling, e*^, be small. Using eqs. (2.2) and 
(2.6) we see that, for finite values of the gauge theory parameters, the inequality e"^ <C 1 can 
only be satisfied up to some critical radius 



5ym 



(2.9) 



Beyond this radius, the M-theory circle opens up to reveal an AdS-j x background (with 
identifications). Similarly, at the corresponding high-energy scales, the five-dimensional 
Yang-Mills theory reveals a UV completion in terms of the (2,0) theory compactified on 
a circle. In the present context, we naturally demand that U„it S> ?7kk, which, using (2.9) 
and (2.6), reduces to 

gt^^-^^l, (2.10) 

where the second inequality follows from eq. (2.8). Eqs. (2.8) and (2.10) imply that the 
supergravity analysis in the D4-soliton background is reliable in precisely the strong-coupling 
regime of the 't Hooft limit of the four-dimensional gauge theory: ^ym — 0, AT^ — oo, S'ym-^c 
fixed and large. 



^This is the tension of a string lying at U — Ukk and extending along x. It can also be computed by 
projecting such a string onto the boundary and integrating (F^) across a transverse section of the string [22]. 
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2.2 D6-brane embeddings 

We are now ready to study the embedding of a D6-brane probe in the D4-sohton geometry. 
Asymptotically (as U oo), the D6-brane is embedded as described by the array (1.1). The 
analysis is greatly simplified by introducing isotropic coordinates in the 56789-directions. 
Towards this end, we first define a new radial coordinate, p, related to U by 



2/3 



and then five coordinates z = {z^ , . . . , z^) such that p = |z| and dz ■ dz = dp^ + p^ d^l^ . In 
terms of these coordinates the metric (2.1) becomes 

/TT\ 3/2 

ds^ = f - J {ri^,u dx^dx" + f{U)dT'^) + K{p) dz-dz, (2.12) 

where 

p3/2r/l/2 

K{p) ^ ^ . (2.13) 

p^ 

Here U is now thought of as a function of p. Finally, to exploit the symmetries of the D6- 
brane embedding we seek, we introduce spherical coordinates A, Jl^ for the ^;^'^'^-space and 
polar coordinates r, 4> for the z^'^-space. The final form of the D4-brane metric is then 



ds^ 



\ R J ^'^^'^ dx^dx" + fiU)dT^) + K{p) {d\^ + A^ d^l + dr'^ + d^^) , (2.14) 



where p"^ = A^ + r^. 

In these coordinates the D6-brane embedding takes a particularly simple form. We use 

x^, A and $7^ (or cr°, a = 0, ... ,6, collectively) as worldvolumc coordinates. The D6-branc's 
position in the 89-plane is specified as r = r(A), (f> = (f>o, where (pQ is a constant. Note that A 
is the only variable on which r is allowed to depend, by translational and rotational symmetry 
in the 0123- and 567-directions, respectively. We also set r = constant in the following, which 
corresponds to a single D6-brane localized in the circle direction. We will study configurations 
in which this condition is relaxed in section 7. 

With this ansatz for the embedding, the induced metric on the D6-brane, gab, takes the 
form 

/f/X 3/2 

ds'^{g) = i^j r]^v dxf'dx" + K{p) [(1 + r^) ^^A^ + A^ d^l] , (2.15) 
where f = dxr. The D6-brane action becomes 

Sd% = -j:^^J d'ae-^^r^g = -T^, j d'a^^h(^ + ^^ \^^l + r\ (2.16) 

where Tog = In j gsi^'nlgf is the six-branc tension and h is the determinant of the metric on 
the round unit two-sphere. Recall that p is a function of A both explicitly and through its 



-7- 



dependence on r(A). Hence the equation of motion for r(A) is 



dX 



1 + 



4p^ 



2 p5 



(2.17) 



Note that r(A) = ro, where ro is a constant, is a solution in the super symmetric hmit {Ukk = 
0), as in [8, 12]. This reflects the BPS nature of the system, which imphes that there is 
no force on the D6-brane regardless of its position in the 89-plane. In particular, then, the 
solution with rg = preserves the U{l)j^ rotational symmetry in the 89-directions. If ?7kk 7^ 
the force on the D6-brane no longer vanishes and causes it to bend as dictated by the equation 
of motion above. We will see below that in this case there are no (physical) solutions that 
preserve the C/(1)a symmetry. 

If C/kk 7^ 0, the analysis is facilitated by rescaling to dimensionless variables as follows: 



A ->■ ?7kkA , r ->■ UjiKT , p ->■ UkkP , 
in terms of which equation (2.17) becomes 



(2.18) 



d_ 
dX 



1 

4/93 



A^ 



We seek solutions such that the asymptotic separation of the D6-brane and the D4-branes, 
L, is finite. That is, as A — 00, r(A) — r^o with 



L 



(2.20) 



In the region A — > 00 we have r 
equation (2.19) with the result 



and p ~ A. Under these conditions we can linearize 



dX 



[A 



^ T*] C± 



3 J. 

2j 



the solution to which is 



_1_ ^ / 2V^3 



-1/2^ 



1 



+ B — Ji/3 



'21/23-1/2^ 
A3/2 



(2.21) 



(2.22) 



where Ji, arc Bcsscl functions and A and B arc arbitrary constants. For large A, the function 
multiplied by A tends to a constant, whereas that multiplied by B behaves as 1/A. This 
means that for large A we have 

r(A) ~ r-oo + I , (2.23) 

where the constants Too and c are related to the quark mass and the chiral condensate, as we 
now show. 
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The bare quark mass is given by the string tension times the asymptotic distance between 
the D4- and the D6-brane.^ Taking (2.20) into account, we have 



27r^2 27r£2 



(2.24) 



The quark condensate can now be computed using a simple argument. The Hamiltonian 
density of the theory can be written as 

n = no + m^J d^OQQ, (2.25) 

where Ho is rriq-independent, and Q, Q represent the hypermultiplet superfields in M = I 
notation. It follows that 

^ = {J d'eQQ) = {^^), (2.26) 

where 6 = (H) is the vacuum energy density. In the last equality, we have assumed that the 
vacuum expectation value of the fundamental scalars vanishes, since they are massive (even 
with TiT-q = after supersymmetry breaking), that is, the energy density increases if they 
acquire a non-zero expectation value. 

Since in the string description the quark mass is given by the asymptotic position of the 
D6-brane, we need to evaluate the change in energy of the D6-brane associated to a change 
of the boundary condition r^o- In view of (2.16) this is given by 

2 



(2.27) 



where C is the D6-brane Lagrangian density rescaled in accordance with (2.18). Using the 
equations of motion this becomes 

\=oo 



S£ = -47r ( dr—— , - ^/.j.d6>^kk \ ^ ^ 773 



df 



A=0 

(2.28) 

where we used (2.23) and the fact that r|A=o = Oj by rotational symmetry in the 567- 
directions. In view of (2.24) we finally obtain 

ii'i^) = J^= -STrH^T^eUl^ c . (2.29) 

As anticipated, the quark condensate is directly related to c; we will see below that c is 
determined by Too- Note that, in terms of gauge theory quantities, 

^(V;V)^5YMiV.M|^c. (2.30) 



^This mass is trivially derived for the brans array (1.1) in asymptotically flat space. With supersymmetry 
{Ukk = 0), this bare mass persists in the decoupling limit and is then inherited by the compactified theory with 

supcrsymmctry-brcaking boundary conditions, since setting Ukk 7^ docs not alter the asymptotic properties 
of the model that, as usual in AdS/CFT-like dualities, determine the gauge theory parameters. 
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The normalization on the left-hand side is that expected in the 't Hooft limit. 

We now return to the full equation of motion (2.19), which we will solve numerically. In 
order to understand the nature of the solutions better, however, it is helpful to first consider 
the limit Too S> 1, in which (approximate) analytical solutions can be found. 

If Too ^ 1, the entire D6-brane should lie far away from the 'bolt' at the center of the 
D4-soliton, hence wc expect the solution to be a small perturbation around the solution for 
Ukk = 0. In other words, we set r(A) = Too + 5r(X) and assume that 6r <C 1. Substituting 
this into (2.19) and expanding to linear order in Sr we find 

_l(A2^r) ~--A2 ^2^^, (2.31) 

which is easily integrated with the result 

<5r~-^^ Tdx ^^^ = --i ^ (2.32) 

Here we have imposed the boundary condition r = at A = 0, as required by regularity of the 
solution at the origin of the 567-space. Integrating once more, with the boundary condition 
'^'^U=oo = 0, we find 

r(A) ~ Too + , ^ . (2.33) 



2roo U2 



2 



From this large-roo solution, we see that the D6-brane bends 'outwards', that is, it is 'repelled' 
by the D4-branes. This repulsion persists for arbitrary values of r^o, as is confirmed by 
numerical analysis. Eq. (2.19) can be integrated numerically for any value of Too-, and we 
have plotted solutions for several values of Too in Figure 1. 

Figure 2 displays the function c{r^) found by numerical integrations. As mentioned 
above, the value of together with the requirement of regularity at A = (i.e., r|A=o = 0) 
determines the solution completely. Hence choosing Too fixes c = c(roo)- Recalling eqs. (2.24) 
and (2.29), which relate r^o and c to the quark mass and condensate, respectively, we see 
that this corresponds precisely to what is expected on field-theoretic grounds: once the quark 
mass is specified, the infrared dynamics determines the chiral condensate. In the field theory 
we expect the U{\)p, symmetry to be spontaneously broken by a non-zero condensate in the 
limit mq 0. This is indeed confirmed by the numerical results for the D6-brane embedding, 
since we see that, in figure 2, c{r oo) approaches a finite constant in the limit r^o 0. In this 
limit, the asymptotic boundary condition on the D6-brane respects the rotational symmetry 
in the 89-plane, but it is energetically favourable for the D6-brane to bend out in the 89- 
plane. Hence the embedding spontaneously breaks this rotational symmetry, which matches 
the C/(1)a symmetry-breaking in the field theory. 

We see from (2.33) that c(roo) ^ l/(2roo) for large Too, which is confirmed by the numer- 
ical results. This implies that the chiral condensate scales as 1/mq for large rriq, as expected 
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Figure 1: The D6-brane embedding for several values of the quark mass oc Too- The interior of 
the 'circle' in the center is the region U < Ukk, which is actually not part of the space. The '2' and 
'3' embeddings correspond to the same value of roo but have opposite-sign c's; the '3' embedding has 
lower energy. 



on field theory grounds [20]. A quick, somewhat heuristic argument for this in QCD is as 
follows.^ The trace of the energy-momentum tensor reads 

T'^^ = m.V'V - ^ TrF^ , (2.34) 

where we have explicitly written down the classical contribution, as well as the anomalous 
one coming from the one-loop beta-function. In the limit rriq — > oo the theory resembles QCD 
with no quarks, for which 

=-11^TVF2. (2.35) 
247r ^ ' 

Taking vacuum expectation values and equating the two expressions we find 
Assuming that (TrF^) 7^ we conclude that (ipip) oc l/niq. 



*We thank D. T. Son for explaining this point to us. 
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For each value of Too > 0, the numerical analysis actually finds two regular solutions, 
which have c's of opposite sign. The curves labelled '2' and '3' in figure 1 are a representative 
example of such a pair of solutions.^ We have confirmed that the solutions with positive c, 
for which r(A) > for all A {e.g., curve 3), have the lower energy density. The negative-c 
solutions bend around the opposite side of the bolt at the center of the D4-soliton geometry. 
Naively, these solutions should be unstable, since the D6-brane is free to 'slide off the bolt' 
by moving out of the (f> = <po plane, and relax to the lower-energy, positive-c solution. We 
will return to this issue in the next section, where we study fluctuations of the D6-brane 
worldvolume fields around embedding solutions of each type — in fact, we will find tachyonic 
fluctuations for the negative-c solutions, confirming the intuition given above. 

0.8 I 1 1 1 1 1 

0.7 - 



C 0.4 




I 1 1 1 1 1 =1 

5 10 15 20 25 30 



Figure 2: The quark condensate {ipip) cx c as a function of the quark mass rriq oc Too- 

Note that r(A) = is an exact mathematical solution of eq. (2.17) even if C/kk 7^ 0. 
Therefore one might be tempted to conclude that it provides an alternative, ?7(l)A-preserving 
solution whose energy should be shown to be higher than that of the J7(l)A-breaking solution 

above in order to establish spontaneous chiral symmetry breaking. However, the solution 
r(A) = is not physically acceptable if U^k 7^ because it corresponds to an open D6- 
brane. That is, this solution yields a D6-brane with a boundary, which is forbidden by charge 
conservation. The boundary is easily seen as follows. With r(A) = we have p = A, so 

^To describe the second solution (curve 2), we allow r(A) < and restrict 4> to the range [0,7r). 
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this solution terminates at the origin of the t/r-plane, that is, at U = [/kk- Since, at this 
point, the in the metric (2.15) still has a finite radius, B?/^Ul/^, the D6-brane would 
have an R^'^ x S'^-boundary if it terminated at [/ = [/kk- Note that if [/kk = then the 
two-sphere radius shrinks to zero-size and the six-brane closes off at C/ = 0. Hence the 
U (l)A-preserving solution r(A) = is physically sensible in the supersymmetric case. This 
result is in agreement with the fact that unbroken supersymmetry forbids a non-zero chiral 
condensate. Note, however, that this is only a formal argument, since in the limit L^kk 
the curvature at U = U^l^ diverges, meaning that the supergravity description cannot be 
trusted in this region [16]. 

In section 7, we will see that the solution r(A) = can play a role in the case [/kk if 
it is simply extended past the origin of the f/r-plane. The physical interpretation of such a 
solution is not that of a D4/D6 intersection, but rather that of a D4/D6/D6 intersection. 

3. Meson spectroscopy from D6-brane fluctuations {Nf — 1) 

Here, we consider a certain class of fluctuations of the D6-brane around the embeddings 
described in the previous section. We consider first the positive-c (negative-condensate) 
embeddings, which were argued to be stable and so should correspond to the true 'vacuum 
state' for a given quark mass. For simplicity we restrict ourselves to fluctuations of the fields 
r and (/>, as these are sufficient to illustrate the physics we wish to exhibit. We are therefore 
considering embeddings of the D6-brane of the form 

= -)- (5^ , r = rvac(A) + 6r , r = constant , (3.1) 

where the fluctuations 5(j) and 5r are functions of all of the worldvolume coordinates and 
rvac(A) is the numerically-determined vacuum embedding of the previous section. 

In corroboration of our previous argument for the stability of these embeddings, we find, 
in particular, that the spectra of these fluctuations arc non-tachyonic. We also briefly consider 
fluctuations around the negative-c (positive-condensate) embeddings and show their spectra 
do contain tachyonic modes, making manifest the instability of these conflgurations. 

In the dual gauge theory, the r- and the (;^-fluctuations correspond to a class of scalar and 
pseudo-scalar mesons, respectively. To see this, consider the complex scalar field parametris- 
ing the 89-plane, X = + iX^ = re^^. Gauge theory Lorentz transformations act only 
on the 0123-directions, under which X is inert, so X-fluctuations clearly correspond to spin- 
zero mesons. However, a parity transformation in the gauge theory corresponds to a ten- 
dimensional transformation that not only reverses the sign of (say) X3, but also replaces X 
by its complex conjugate, X] this leaves r invariant but reverses the sign of 0. Prom the gauge 
theory viewpoint, this is clear from the presence in the microscopic Lagrangian of a coupling 
of the form tp^Xi/jj^ + t/jlXipT^: the transformation X3 —>■ —X3 exchanges V'r and V'l, so this 
term is only invariant if, simultaneously, X ^ X. Prom the string theory viewpoint, the fact 
that only the combination of both transformations is a symmetry can be seen by examining 
the fermionic spectrum of the D4-D6 open strings. This arises from the zero modes of the 
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world-sheet fermions in the Ramond sector, that is, in the NN and DD directions: the 23- 
and 89-directions (in the hght-cone gauge) in (1.1). Each of these modes is labelled by its 
weight under rotations in each of the 23- and 89-planes, (si,S2),Sj = ±1/2. The GSO pro- 
jection requires si = —S2- Since under X3 — > — X3 we have (si,S2) (— si,S2), and under 
Xg —Xg we have (si,S2) (si,— S2), we see that the spectrum is only invariant under 
the combination of both transformations. 

3.1 Analysis of the spectrum 

Using the background metric (2.14) with the rescahngs (2.18), the pullback of (2.14) to the 
embedding given in (3.1) is 

3/2 

= / - j Tfi^^dx^dx'' + K[{1 + r^JdX^ + X^dnl + 2n,,5„((5r)dAdx"] 

+ K[daiSr)db{5r)dx''dx'' + (r,,, + 5rfda{d(j))db{6<p)dx''dx''] , (3.2) 

where a and b run over all of the worldvolume directions. To quadratic order in fluctuations, 
the D6-brane Lagrangian density is then 



L-A^)^3(4.L-A^) ^^^^ 



16pio 4p7 



3rvac f-^^ 1 ^ r.^dx{5r^) 



Ap^ \ 4p3 7 l + f 

rva,c \ r vac / ' ■ 

2 /\Pl (K^W2 



where £0 is the Lagrangian density evaluated for the vacuum embedding, as in eq. (2.16). 
Here, = -|- and the gaa in the last line correspond to the metric coefficients from 
eq. (2.15). In particular then, they contain no dependence on the fluctuations. Of course, 
integration by parts and the equation of motion (2.17) for r^ac allowed the terms linear in Sr 
to be eliminated. 

The linearised equations of motion that follow are, for b(^: 

a J.2 / 1 \4/3 „2 / 1 \2 

= 0, (3.4) 



4M2 



1 d 



1 + dxm 



+ ''vac V 4^vac 



and, for 5r: 
A2 



4M2 o' 



vac \ vac / \ r vac / 
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(1 + ^0)=^/' V ^p- 



vac 



2 



d_ 
dX 



* ' vac ' vac / -j^ -■- 



,3 

vac / 



5r 



In these expressions, is the Laplacian on the two-sphere. 

Prom the form of the equations of motion, it is clear that we may separate variables to 
write 

5cl> = V{\)e'^'^-^ Y,^m,{S'') , Sr = 7^(A)e^'='•- Yi^mAS') , (3.6) 

where are the conventional spherical harmonics. We then look for normalizable solutions 
with definite S'^-angular momentum i^^^ and four-dimensional mass M| ^ = —k'^j.- For the 
sake of simplicity in the following, we only focus on the case ^^^j. = and find the lowest-mass 
modes for each fluctuation. 

We employ the shooting method to solve the linearized equations of motion. First we 
use the asymptotic form of the equations for A 00 and the requirement of normalizability 
to determine the asymptotic boundary conditions for the fluctuations. Then we vary the 
four-dimensional mass parameters until, by numerically integrating the equation of motion, 
we find solutions that are also regular at the origin. Recall^ that the asymptotic behaviour 
of the vacuum profile is r^^^ ~ r^o + c/A, where r^o is related to the quark mass by equation 
(2.24). In this asymptotic region we have that Pvac(A) A, so the equation of motion for S(f) 
becomes, approximately, 

9M? r2 r] 

and similarly that for 5r becomes 

^^,^,) + ^[A^a,M=0. (3.8, 

It follows that the asymptotic behaviour of 5r(X) docs not depend on that of rvac(A), whereas 
that of d(p{X) depends crucially on whether or not Too = 0, that is, whether or not the quarks 
are massless. We will now show that if = then there exists a normalizable Scj) mode with 

= 0, whereas if Too 7^ there is a mass gap in the spectrum. 

If we set Too = 0, we find that the asymptotic equation for becomes 

1 

:(6ct>)+dlid<P)=0. (3.9) 



4M2^ A3 



The first term is sub-leading for large A, so we may drop it to find that the asymptotic form 
of 6(f), regardless of the value of M^, is 

6(t)c^b + a\, (3.10) 



^Also recall that we are working with dimensionless variables r and A, rescaled as in eq. (2.18). 
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Figure 3: Squared masses for the four lowest-lying S(p modes (dashed lines), dual to pseudo-scalar 
mesons, and the four lowest-lying dr modes (solid lines), dual to scalar mesons, as functions of the 
quark mass niq oc Too- 



where a and b are arbitrary constants. For r^o 7^ 0, however, we find that 

t2 



9M^ 1 2 



4M2 A3' 



A 



0. 



Again the first term is sub-leading for large A, so the solution is 

S(f) a + b/X . 



(3.11) 



(3.12) 



Normalizable solutions correspond to the boundary conditions a = and a = in the cases 
''oo = and Too 7^ 0, respectively, since the terms with coefficients b and b are the ones that 
fall off more rapidly at infinity. Since the differential equations we are solving are linear, we 
may choose 5 = 6=1 without loss of generality, so we are left with the boundary conditions 
{d(P = 1 , 6(j)' = OjU^oo if Too = 0, and {(50 = 1/A , 6<j)' = -l/A^jU^^ if 0. A 

similar analysis for the Sr fluctuations shows that they obey the same asymptotic equation 
of motion as the S(f) fluctuations in the Too 7^ case, so the appropriate boundary conditions 
are {^r = 1/A , 5r' = -l/A^jU^^. 

Now note that Sep = e**^'*'^, with k"^ = 0, is an exact, regular solution of (3.4) with zero 
four-dimensional mass (and V^{6(p) = 0), which is normalisable only if Too = 0. This means 
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Figure 4: (a) Linear fit of the lowest-lying 5(j) mode for different values of rriq ~ Too- (b) Degeneracy 
of the lowest-lying Scj) (dashed) and Sr (solid) modes in the supersymmetric limit Too 1- 

that with masslcss quarks there is a massless, pseudo-scalar meson in the spectrum of the 
gauge theory, whose dual is the zero-mode fluctuation of the D6-brane field (p. The math- 
ematical form of these fluctuations shows that they correspond to rotations of the vacuum 
embedding in the 89-plane. Of course, this form is precisely as expected for the Goldstone 
mode associated with the spontaneous breaking of this rotational symmetry. 

The shooting technique allows us to verify this numerically, and to find the values of 
M^r for which other normalizable, regular solutions exist. The results are the low-lying 
(pseudo)scalar meson spectra, which arc displayed in Figure 3 for different values of the 
quark mass. As anticipated, there is no mass gap in the (?!)-meson spectrum when the quark 
is massless. There is a smooth transition to a theory with a mass gap as the quark mass 
is increased and, in fact, a best-fit (shown in figure 4a) yields a linear relationship,''' M| ~ 
0.73 Too, for small expected from the GMOR relationship (1.2), which we 

will reproduce analytically in the next section. 

As shown in figure 3, there are no normalisable, regular solutions of equation (3.5) with 
Mr = 0, so there is a mass gap in the r-meson spectrum regardless of the value of the quark 
mass. For small quark mass the mass scale of all of the mesons, except the lowest-lying 5(j) 
mode, is^ 

M'-M^r.^. (3.13) 

This is another reflection of the lack of decoupling between the QCD scale and the compact- 
iflcation scale. 

'^An analogous rclationsliip was observed in [14]. 

®Tliis also sets the scale of the glueball mass spectrum [3]. 



-17- 



As suggested by Figure 3, the r- and ^-mesons occur in pairs that become degenerate for 
large enough a quark mass. The reason is that, in the supersymmetric case (C/kk = 0); these 
mesons belong to the same supermultiplet and hence have equal masses, as in [12]. If C/kk 7^ 
this degeneracy is only approximate for quark masses much larger than the supersymmetry- 
breaking scale, M^nsy, becoming exact in the limit niq oo. Geometrically, this corresponds 
to the fact that the distortion of the D6-brane away from the flat, supersymmetric profile 
decreases the farther away from the D4-branes we embed it. Since this distance scale is 
measured by L oc m^, and C/kk characterises the supersymmetry breaking, we expect that, 
for L S> C/kk, the characteristic mass scale of the mesons should become independent of C/kk 
and depend, instead, on some ratio of L and R. To see that this is indeed the case, we 
approximate the vacuum profile by the flat solution in equations (3.4) and (3.5), that is, we 
set rvac(A) = Too, and we expand these equations assuming that Too 3> 1. In terms of a new 
variable y = \/L we obtain (for zero S'^-angular momentum) 



L(2/2 + 1)3/2 ^ ' y^dy 

and 



+ ^:rbW] = o> (3-14) 



;Sr + ^-^[y%6r]=0. (3.15) 



L(y2 + 1)3/2 y2dy^ 

The fact that the two types of fluctuations satisfy the same equation confirms that the spec- 
trum becomes degenerate in this limit. It is also clear that the mass scale for the mesons in 
this limit is now 

since this is the only scale that appears in eq. (3.14), (3.15). 

The transition between the regimes (3.13) and (3.16) occurs at L ~ C/kk, and is illustrated 

for the lightest r- and 0-mesons in figure 4b, which was generated by solving equations (3.4) 
and (3.5) numerically. The figure clearly illustrates the degeneracy of the modes in the 
expected regime. The condition L ~ C^kk translates into rriq M^,,^y, with M^^^y as defined 
in (1.4), so Mg^gy is indeed the scale at which supersymmetry is restored, as suggested by 
our choice of notation. Note that at the transition point, at which of course eqs. (3.13) and 
(3.16) coincide, we have M ~ m^l{g'^^N^). As we are working in the 't Hooft limit of the 
gauge theory, these mesons are very deeply bound, i.e., M <^ 2m^. The suppression of these 
meson masses by the 't Hooft coupling is even stronger than that found in a similar context 
in [12], where M ~ wiq/VwYM^- -^^^ large quark mass there is an extra suppression implicit 
in eq. (3.16) due to the fact that scales linearly in m^, as opposed to quadratically. 

As we are working in a range where we expect the gauge theory is five-dimensional, it is 
perhaps more appropriate to re-express the transition scale in terms of the five-dimensional 
gauge coupling; it then reads ~ g'^N^M^^. Similarly, the mass scale typical of the mesons 
becomes ~ m^/ {g'j^N^). In any event, the transition scale suggests that the effective 
mass squared of some of the microscopic degrees of freedom must be at least of the order 
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M^nsy ^ -^KK, father than just Mkk- Note that this mass scale is a factor of V'^ym-^c larger 
than that coming from a simple one-loop calculation for the scalar masses, which yields 
V^YM-^c-^KK- This factor seems to be ubiquitous in extrapolations from weak to strong 't 
Hooft coupling in the context of AdS/CFT. 

As we saw above, the squared masses of all mesons scale linearly in the quark mass as 
rriq — >■ oo. For the lightest meson, the results displayed in figure 4b seem to indicate that, in 
fact, this linear relationship is valid not only in the limits of small and large rriq, but for all 
quark masses. This result has been confirmed to within the best accuracy of our numerical 
analysis. Unfortunately, eq. (3.14) has four regular singular points and so no closed-form 
analytic solution is available, but its numerical solution gives, for the lowest-mass modes, 

~ 1.66 ~ 20.91 4%-"^<J - 0-74 M^^^oo • (3.17) 

The final form of this expression agrees with a best-fit of the data plotted in the region of 
figure 4b with ^ 1 and, in principle, describes the extension of this figure to arbitrarily 
high values of r^o- We therefore conclude that, in our model, the GMOR linear relation 
between and extrapolates to all quark masses. 

3.2 Stability Issues 

Recall that, at the end of section 2, we remarked that our numerical integrations had revealed 
two classes of six-brane embeddings, those with positive c and those with negative c. From the 
numerical analysis, we could show that, for a given Too, the positive-c or negative-condensate 
solutions had the lower energy density. Further we argued that the negative-c solutions should 
be unstable. 

The stability of the negative-condensate (positive-c) vacuum embedding is reflected in 
the fact that the spectrum of fluctuations does not contain tachyonic modes. While numerical 
searches did not reveal any tachyonic fluctuations, one can formulate a general argument that 
no such modes can exist in the spectrum of 5(j).^ Consider the equation of motion (3.4) (with 
zero angular momentum on the two-sphere, for simplicity). Multiplying by d(p, integrating 
over A and integrating by parts, yields 




Note that we were able to set the term on the right-hand side to zero using the behaviour of 
the normalizable Scf) modes at A = and A = 00. On the left-hand side, the second term in 

®We expect that these arguments extend to Sr fluctuations but we did not expUcitly consider these modes 
a possible decay channel between the two D6-brane embeddings, due to the geometric obstruction provided 
by the 'bolt' at the centre of the D4-soliton. 
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the integrand is manifestly negative or zero, while the sign of the first term depends on that 
of M|. Therefore, for tachyonic modes it is clear that the integral will be negative-definite, 
which is inconsistent with the vanishing of the right-hand side. The argument extends to 
include nonvanishing which simply introduces an additional negative-definite term under 
the integral. Hence we conclude that the negative-condensate embcddings are free from 
tachyonic instabilities in this sector. As this was the sector with the lowest lying mode, we 
might have expected it to be the most potentially problematic with regards to stability. Hence 
we are confident that these positive-c solutions are stable. 

Now we turn our attention to the negative-c or positive-condensate embeddings. Recall 
that our intuition was that these solutions should be unstable to 'sliding off the bolt.' Since 
the embedding coordinate (r, here) always vanishes at some value of A for this case, we 
cannot address the issue of such an instability using the dcp equation of motion, since (p is 
not everywhere well-defined. Hence, we change to Cartesian coordinates X and Y in the 
r(?!)-plane. With the choice that the D6-brane lies at <;livac(A) = (f)o = {i.e., l^ac(A) = 0), it 
follows that the embedding profile in X satisfies the equation of motion (2.19) and we may 
set X^^^{X) = r,.ac(A). Fluctuations in Y{X) are now the modes relevant for stability and have 
the benefit of being well-defined for all A. One finds the following equation of motion for 
these fluctuations: 



1 



dX 



A2 / 1 ^2 



y2 

vac 



+ ll±i^i)iy = 0, (3^9) 



where, now, p^^^ = A^ + X"^^^. Making a stability argument along the lines of that given 
above for 5(f) is no longer possible, due to the presence of the final term in (3.19). This term 
introduces into the integrand analogous to that in eq. (3.18) a contribution that is manifestly 
positive. Hence this approach does not yield a definite conclusion in this case. 

However, it is still possible to look directly for tachyonic modes in the spectrum of fluc- 
tuations around the negative-c embeddings using our numerical techniques. Prom eq. (3.19) 
one can deduce the asymptotic behaviour of the 5Y modes to be 

(5y ~ a + 6/A , (3.20) 

independent of the vacuum embedding. Using the shooting method as before, the spectrum 
of each positive-condensate embedding considered numerically is found to contain a tachyonic 
mode and, furthermore, these modes become 'more tachyonic' as the quark mass increases — 
see figure 5. As is suggested in the figure, the squared mass of the lowest-lying SY mode has a 
linear dependence on the quark mass. In fact, a best-fit of the data gives ~ -0.72 M| 



An examination of the next-to-lowcst-lying 6Y mode (not shown in the figure) indicates 
that the mass of this mode, which is real and positive for rriq = 0, also decreases as a function 
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of Too- So for sufficiently large quark mass the spectrum may contain more than one tachyon. 
However, we could not confirm this directly because for large the negative-c embeddings 
approach extremely close to the bolt, and our computer code behaved erratically in this 
situation. 



4. The pseudo-Goldstone boson 

As discussed in the previous section, if rriq = 
the spectrum contains a massless, pseudo- 
scalar mode that can be understood as the 
Goldstone boson of the spontaneously bro- 
ken U{1)a symmetry. This is akin to the 
rj' meson in QCD, which is massless in the 
large- limit. In an abuse of language, how- 
ever, we will refer to this Goldstone boson 
as a 'pion'. For non-zero quark mass, the 
U{1)a symmetry is explicitly broken and the 
pion becomes a pseudo-Goldstone boson with 
a mass M^- that, according to our numerical 
analysis, grows as ~ m^. This means 
that at low energies, E <^ Aqcd, the dynam- 
ics is dominated by this particle. In this sec- 
tion we will provide an analytic proof, using 
the string description, that the pion's mass 
actually obeys the well-known Gell-Mann-Oakes 
Renner relation (1.2). 

Recall that the gauge theory pion is dual to a 




Figure 5: The lowest- lying dY mode for some 
negative-condensate embeddings (upper line), and 

the lowest-lying SY mode for the corresponding 
positive-condensate embeddings (lower line) — 
the SY modes are tachyonic. 



-fluctuation of the D6-brane of the form 



Sep = e </?(A) that solves equation (3.4) with V {6(j)) = 0, that is, 

\p{X)^] = -MV(A)¥' , 
where = —k'^ is the four-dimensional mass. 



(4.1) 



Ai(A) = X' 



\2 2 

A r 

vac 



1 + 



4p: 



ViT^ 



4:ML 



1 + 



4p: 



4/3 



and pI^^ = a -|- r^^^. The allowed solutions are those that both satisfy (^U=o 
normalizable with respect to the scalar product 



(4.2) 



(4.3) 



and are 



Jo 



dA/x(A)<^t(A)<^2(A). 



(4.4) 
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This defines an eigenvalue problem whose solution is a complete set of orthogonal eigenfunc- 
tions with eigenvalues M^. Notice that if (^(A) = constant then the eigenvalue equation 
is satisfied with = 0, but the mode is normalizable only if r^^^{\) — as A — oo. If 
(^(A) is not constant, multiplying eq.(4.1) by 99, integrating both sides with respect to A and 
performing an integral by parts in the left-hand side, we arrive at 

2 rf°dAp(A)<^2 

/„°°dAM(AV2>"- ^^-'^ 

Since from (4.2) and (4.3) we have ^»(A) > and > 0, it follows that is positive, as 
indicated. ^° 

Having said that, we now want to compute the mode with lowest eigenvalue > when 
mq is non- vanishing but small, that is, let us suppose that r^ac — ^ foo A — > 00, with Too 
small but non-zero. In this case we expect that the normalizable cigcnmodc can be obtained 
from the massless-quark zero- mode, namely V'(A) = constant, by using perturbation theory. 
To this end, it is convenient to define 



^(A) = Vp(A)¥'(A), (4.6) 
which satisfies the new eigenvalue equation 

^_|^ = _m2i.V, (4.7) 

where ^ = ^/p, and u = ^/pis the new integration measure. If r^o = then ipiy^) = ^(A) is a 
normalizable solution of (4.7) with zero eigenvalue. Making a small change in the boundary 
condition. Too Sroo ^ 0, induces a small change in the vacuum profile rvac(A), which in 
turn induces a small change in the functions *(A) and ^'(A) appearing in equation (4.7). 
Standard quantTim mechanics perturbation theory can then be applied to obtain the new 
lowest eigenvalue. 

Let us denote with a 'bar' the quantities corresponding to Too = 0, so that for ^ we 
have 

* = * (5^r , (4.8) 
u = u + Su, (4.9) 

where 6"^ and diy arc the differences induced by the change in rvac(A). The new lowest 
eigenvalue and its associated wave-function, > and ijj = + Stp, obey, to leading order, 

• I' f ^\ - - 

5'il;^-^Stlj-Sl-\qi = -M^i?^. (4.10) 



^"This is not valid when rvac vanishes at some point (which is not the case studied here) — see the more 
detailed discussion after equation (3.18). 



- 22 - 



Decomposing 6ijj as a linear combination of massless-quark eigenfunctions, dip = J2^=o 
we obtain 

oo / T \ 

^a,M2-^, + 5 -U = m2/?*, (4.11) 



n=l 



where we have dropped the first term in the series, using Mq = 0- Multiplying by inte- 
grating over A and using the orthogonality condition dA i>'^tpn = for n > (since ^ is 
the zero-mode eigenfunction and is real) we deduce that 



M 



r°° -o r^x -r, r°° _ .. ^ r_ . ^ -ia=oo 

/ d\v^^ = d\5\-\^!^=l d\{^5^-^5^)= ^5^-^5^ -(4.12) 
Recalling that 

^ = Af..__ ( 1 + 

4p: 



* = Af... ( 1 + ^ ) (1 + flj-y^ , (4.13) 



which implies 

5^ = \5r.^ (l + ^) (1 + fL)-V4 _ ^(1 + flj-y^ 5r.^ - 

-Af_(^l + ^) ^{l + flj-^/^f^^5f^^, (4.14) 

we see that the only contribution to the last term in equation (4.12) comes from A = oo. This 
is easily evaluated, since for large A we have 

f.,, -^ + , * ~ Af.,, ~ c + O . (4.15) 

It follows that the new eigenvalue is given by 

,,2 cSvoo c5r^ 

To show that this formula is precisely the GMOR relation (1.2), we need to identify the pion 
decay constant in the string description. In this description, a shift of the chiral angle that 
parametrizes the space of vacua related to each other by the C/(1)a symmetry corresponds 
to a rigid rotation of the D6-brane field (f) (both take values between and 27r). Similarly, 
the pion field is identified with the normalizable, massless mode 54>, described in section 3.1. 
The pion decay constant can be read off from the normalization of the kinetic term in the 
four-dimensional low-energy effective Lagrangian for this mode. That is, we integrate the 
5(j)-\,e]:ms in eq. (3.3) over A and the coordinates on 5^ and as a result find a four-dimensional 
action of the form 

S = 1 d''xd^{5<t>)d^'{5ct>) , (4.17) 
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where /jr is the pion decay constant. From eq. (3.3), then wc have 

poo 

= 47rT^eUl^ / dA/2(A) . (4.18) 
Jo 

Combining this result with eqs. (2.24) and (2.29), we see that eq. (4.16) is precisely the GMOR 
relation (1.2), as anticipated. 

We close this section by verifying the agreement between the analytical expression (4.16) 
above and the results obtained purely by numerical methods and displayed in figure 4a. The 
integral in equation (4.16) can be evaluated numerically with the result 

I ''-^''^-TwJ. ''^{'^w^ ^--w. 

Similarly, the numerical result for the constant c for zero quark mass (see Figure 2) is c ~ 
0.71. Substituting these two values in equation (4.16) we get ~ 0.73 M^j^^roo, in exact 
agreement with a best-fit of the points in figure 4a. 

Recall from section 2 that for nonvanishing quark mass we found that there are two 
possible values of the quark condensate, differing in sign, but that only the negative condensate 
is stable. In section 3.2 the positive-condensate solution was explicitly shown to be unstable 
by showing that the pseudo-Goldstone mode 5Y was tachyonic for this configuration. Since 
the field-theoretic arguments establishing the GMOR relation apply regardless of the sign of 
the quark condensate, the squared mass of this mode should also obey the GMOR relation, 
at least approximately, with an appropriate sign change. Recall that eq. (1.2) represents the 
first term in an expansion around = 0. Now in the string description at mq = 0, the 
D6-brane configurations for positive and negative c are identical, except for a rotation of cf) 
by TT. Hence up to a crucial sign, the numerical coefficients in eq. (4.16) are unchanged and 
hence the tachyonic mode should satisfy My — — 0.73 M^j^fJroo- From Figure 5, it does indeed 
appear that the slopes of the lowest-lying dcj) and 5Y modes are equal, but opposite in sign. 
More precisely, our best-fit of the numerical results for negative c gave My — —0.72 M^^r^, 
which gives agreement to within about 1.4%. 

5. Multiple flavours (Nf > 1) and a holographic Vafa-Witten theorem 

Consider a vector-like gauge theory with vanishing 6'-angle (like QCD) and N( flavours, all with 
identical mass > 0. The Vafa-Witten theorem states that the vector-like U{Nf) flavour 
symmetry cannot be spontaneously broken [19]. This theorem does not exclude the possibility 
that, at mq = 0, the ?7(iVf )-preserving vacuum becomes degenerate with one or more U{Nf)- 
breaking vacua. In this section we discuss the holographic realization of the Vafa-Witten 
theorem, within the approximations implied by the validity of the supergravity/Born-Infeld 
theory. In particular, we are considering iV^ — > oo and A^f <C A^c- Note that, strictly speaking. 
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the Vafa-Witten theorem need not apply to the gauge theory on the D4/D6-brane system 
because of the presence of Yukawa-Uke couphngs between the fundamental fermions and the 
adjoint scalars. However, the latter scalar fields are massive and the theorem should certainly 
apply in a limit where the scalar masses are sent to infinity. Further, one expects the same 
result for large but finite masses since the vacuum should still not be sensitive to the scalars. 
In the present case, the scalar masses are likely of order Mkk ~ ^qcd a-^d so these couplings 
would not be suppressed. Hence, from the gauge theory viewpoint, it is not clear a priori 
whether the effect of these fields can change the vacuum structure to the extent of violating 
the Vafa-Witten theorem.^^ The results of this section answer this question in the negative. 

The U{N() global symmetry of the gauge theory corresponds, in the string description, 
to the U (Nf) gauge symmetry on the worldvolumc of A^f D6-brane probes in the geometry 
(2.1). Their dynamics is described by the so-called non-Abelian Dirac-Born-Infeld (DBI) 
action [23]. This is a generalisation of the single-brane Abelian theory (2.16), in which both 
the gauge fields and the transverse scalars are promoted to Hermitian, U (iVf) matrices. This 
means, in particular, that the positions of the D6-branes are no longer commuting quantities in 
general. In the particular cases in which the transverse scalar matrices can be simultaneously 
diagonalised, each of their A^f eigenvalues can be interpreted as the position of one of the 
D6-branes along the corresponding transverse direction. 

For the U(Ni) global symmetry to be realized in the gauge theory, the mass matrix for 
the fundamental quarks must be diagonal with all eigenvalues equal to m^. In our string 
picture, this translates into the boundary condition that, asymptotically, all the D6-branes 
should have the same asymptotic boundary conditions. That is, they are all aligned with 
each other in the (z^, . . . ,z^) subspace and all lie at a distance 27r£^mq from the origin. We 
would like to show that the configuration (which obviously satisfies the boundary condition 
above) in which all the D6-branes are exactly coincident, and are embedded exactly as in the 
case of a single D6-brane, is the minimum-energy solution of the equations of motion that 
follow from the non-Abelian DBI action. 

To avoid problems with the coordinate singularity at r = 0, we adopt the Cartesian 
coordinates X and Y introduced in section 3.2 to describe the transverse positions of the 
D6-branes in the r^-plane. Then the configuration for the Nf coincident six-branes may be 
written as: 

^(A) = r.ac(A) • I , Y{X) = 0, r = To-I, (5.1) 

where rvac(A) is the profile for a single D6-brane determined in section 2 and I is the Nf x 
Nf identity matrix. Since all the transverse scalars commute, there is no ambiguity in the 
interpretation of this solution as Nf overlapping D6-branes lying along the curve ^(A) = 

rvac(A). 

Let us first argue that this configuration is a stable solution of the non-Abelian DBI 
equations of motion. We start by expanding the non-Abelian DBI action to quadratic order 
in fluctuations around the above configuration. Since there is an overall single trace in front 

^^We thank A. Nelson for a discussion on this point. 
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of the action, to linear order only variations proportional to the identity generator couple to 
the U{1) background fields, as given above. Hence the resulting expression for these fields is 
exactly that obtained from the linear variation of the action for a single D6-brane. It thus 
follows that the configuration (5.1) is a solution of the non-Abelian equations provided rvac(A) 
solves the corresponding Abelian equations. 

We further observe that, to quadratic order in fluctuations, the non-Abelian DBI action 
reduces to the sum of identical copies of the Abelian action, one for each of the 
independent fluctuations. The reasons for this are two-fold: First, given that the background 
flelds in (5.1) are proportional to the identity, the quadratic fluctuations must again contribute 
a term in the overall U{1). Second, the difference between the non-Abelian and the Abelian 
DBI actions consists entirely of terms that involve commutators of fields. Hence, with the 
given background, the only possible nonvanishing commutators will involve two fluctuations, 
but such expressions are necessarily SU (A''f)-valued, i.e., traceless. It follows that upon tracing 
over the gauge indices (and integrating by parts), the quadratic action will take the form 



where Z'^ = X,Y, a is the gauge index and Oij is the same operator as appears in the 
quadratic expansion of the Abelian action. Thus the linearized equations of motion, and 
hence the spectrum, arc just copies of those in the Abelian theory describing a single D6- 
brane embedding. It follows that the solution (5.1) is stable provided its Abelian analogue is 
stable, which we know to be the case from the analysis of section 3. 

Note that these fluctuations of the D6-branes, which give rise to the copies of 
the spectrum, are distinct physical excitations and not gauge-equivalent: each of them is 
associated to a different generator of U{Nf), so there is no small U{Nf) gauge transformation 
(that is, one that tends to the identity at infinity in the 0123-directions) that takes one mode 
into another. In particular, this implies that in the limit there are massless, 

pseudo-scalar particles! We will come back to this point in the last section. 

The arguments presented so far show that (5.1) provides a stable extremum of the non- 
Abelian DBI action, and therefore a local minimum of the energy of the Nf D6-branes. Strictly 
speaking, this does not exclude the possibility that there exist other minima. If this was the 
case, then the true vacuum (or vacua) would be that (or those) with the lowest energy. On 
physical grounds, however, the existence of these other vacua is unlikely. The reason is as fol- 
lows. If we ignore the interactions of the D6-branes with each other, then it is clear that (5.1) 
is the minimum-energy solution amongst those satisfying the boundary conditions described 
above, since in the absence of interactions each D6-brane must minimise its own energy. If 
we include open string tree-level interactions by describing the D6-branes' dynamics with 
the non-Abelian DBI action, then the only obvious configuration that satisfies the appropri- 
ate boundary conditions is (5.1), which we have shown to be a stable solution. Inclusion 
of one-loop (or higher) open string interactions between the D6-branes (such as exchange of 




a=l 



(5.2) 
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closed strings) would presumably just strengthen the conclusion that (5.1) is the preferred 
solution. To see this, consider the simple case of TVf = 2 D6-branes that overlap with each 
other asymptotically, but that are slightly misaligned in some small region. Locally, this 
situation is analogous to that of two planar D6-branes in flat space that arc almost but not 
exactly parallel, in which we know that there is a net force that tends to align the brancs. It 
seems plausible that a similar force will act on the brancs in the situation of interest here, and 
therefore that the preferred configuration will be (5.1), in which all the branes are perfectly 
aligned with each other. 

This conclusion provides a holographic realisation of the Vafa-Witten theorem, since 
the U{Nf) gauge symmetry on the D6-branes is unbroken if and only if they are all exactly 
coincident. 

By continuity, the arguments above imply that, in the limit 0, the J7(A''f)-preserving 

configuration has an energy no greater than that of any other configuration. Just like in the 
Vafa-Witten theorem, however, this does not rule out the possibility that in this limit some 
other state becomes degenerate with that in eq. (5.1). In fact, the appearance of massless 
modes in the linearized spectrum at = suggests that the D6-branes might be free to 
rotate independently in the 89-plane.^^ To see that this is indeed the case, consider, for 
simplicity, the N( = 2 configuration 

X{X) = rvac(A) cos 00 • I , ^(A) = r^e^{X) sin 0o • 0-5 , r = tq • I , (5.3) 

where as = diag(l, —1) and (po is some fixed angle. The radial vacuum profile above is that 
corresponding to Too = 0. This configuration consists of two D6-branes with the same radial 
profiles but separated an angle 2(f)o in the 89-plane, and obviously satisfies the boundary 
condition that the branes overlap at A — > 00. To see that it solves the non-Abelain equations 
of motion, consider expanding the non-Abclian DBI action to linear order in fluctuations 
around (5.3), as before. Since commutator terms are irrelevant at this order, and X and 
Y are diagonal, the non-Abelian action again reduces to the sum of two identical Abelian 
actions. In fact, the same argument shows that not only is the configuration (5.3) a solution 
of the non-Abelian equations, but also that its energy is, for all (f)o, equal to twice that of 
a single D6-brane. It follows that (5.3) is stable for all 0o. To see this, suppose that it is 
not, i.e., that there is some fluctuation around it that is tachyonic. Presumably, this would 
mean that this solution can decay to another one of lower energy. However, the latter would 
also have lower energy than the solution with cpo = 0, in contradiction with our continuity 
argument above. 

Note that the previous reasoning applies equally well to any distribution of Nf D6-branes 
around the (^circle, since it only relies on X and Y being diagonal. Thus, at = 0, 
all these configurations are degenerate, lowest-energy, stable solutions. Although it may be 
an interesting exercise to verify their stability explicitly by expanding the non-Abelian DBI 

^■^Notc that this is not guaranteed by the existence of massless modes, since this does not exclude the 
appearance of higher order potential terms for the fluctuations, e.g. . 
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action to quadratic order in fluctuations around (5.3), one should keep in mind that this 
conclusion is only expected to hold at leading order in the ~ Qg expansion, since, as 

discussed above, we expect higher-order effects to lead to an attractive force between the 
branes, which would tend to align them. 



6. Finite temperature physics 

The gauge theory we have so far been describing is at zero temperature. To study the theory 
at a finite temperature T, we may follow the standard prescription: analytically continue the 
time coordinate t ^ = it, periodically identify with period dt^ = 1/T, and impose anti- 
periodic boundary conditions on the fermions around the tE-circle. This prescription applies 
equally well in the gauge theory and in the dual string description. In the latter, however, 
one must consider all possible spacetime solutions with the appropriate boundary conditions, 
that is, possible saddle points of the Euclidean path integral over supergravity (or rather, 
string) configurations [2]. 

In the present case, there are two^^ Euclidean supergravity solutions which must be 
considered. The first is simply the Euclideanised version of that appearing in eq. (2.1), 

which dominates at low temperatures. Recall that / = 1 — U^v^/U^- The second is the 
Euclidean black hole coming from a nonextremal D4-brane throat 



ds' = 



which dominates at high temperatures. In this case, / = 1 — U^/U^. In the Lorentzian- 
signature solution, the event horizon is located at U = 11^- 

Of course, these Euclidean metrics are identical up to interchanging the role of and r 
(and replacing C/kk with fZ-r)- To match the boundary conditions set by the finite-temperature 
gauge theory, we identify these coordinates with periods 

in both solutions. Since we wish to avoid conical singularities in either spacetime at U = 
Ukkj Ut, the boundary conditions fix the metric parameters as 



^^In analogy with [24], there may be an infinite family of Euclidean solutions labelled by two integers 

specifying which cycle of the fET-tonis shrinks to zero-size in the bulk. However, given that the ts- and T-axes 
are orthogonal, here either the solutions in eq. (6.1) or (6.2) will dominate the path integral. 
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As we commented above, the metric (6.1) dominates the path integral at low tempera- 
tures, whereas that given by eq. (6.2) dominates at high temperatures. One determines the 
transition point by comparing the Euclidean action of these two solutions [2]. As the met- 
rics describe the same geometry, it is easy to see that the transition occurs precisely when 
St = Sts- That is, a phase transition occurs at the critical temperature 



Prom the gauge theory viewpoint, this is a confinement/deconfinement phase transition. This 
can be seen by the fact that the temporal Wilson loop (around tg) vanishes for the low- 
tcmpcrature background (6.1), whereas it is nonvanishing for the Euclidean black hole (6.2) 
[2]. Alternatively, on the Lorcntzian sections, the Wilson loop as computed from the solution 
(2.1) exhibits an area law behaviour, while that computed for the Lorcntzian black hole does 
not because the string can break in two pieces, each having an endpoint at the horizon [26]. 
One might also observe that, in the low-temperature phase, one has a discrete spectrum of 
glueballs [3], but this gives way to a continuum of excitations in the high-temperature phase 
dual to the black hole background. The analogous phase transition in 3+1 dimensions was 
discussed in ref. [25]. Finally, note that the phase transition takes place at a temperature, 
(6.5), where the gauge theory is starting to behave five-dimensionally, rather than four- 
dimensionally. 

The preceding discussion refers only to the physics of the (compactified) five-dimensional 
gauge theory and is independent of the presence of fundamental matter fields. In previous 
sections we have exhaustively analysed the physics of dynamical quarks in the confining 
phase. We can extend the discussion to the deconfining phase by considering probe D6- 
branes in the black hole background (6.2). In reference [14] the study of D7-branes in an 
asymptotically AdS^ black hole background revealed the existence of two sets of regular 
embeddings, parametrised by the quark mass and chiral condensate of the dual gauge theory, 
and distinguished by whether or not the D7-brane intersects the horizon. It was suggested 
that the gauge theory undergoes a phase transition, signalled by a discontinuity in dc/dm^, 
when the D7-brane undergoes its 'geometric' transition. Here we consider similar physics in 
our four-dimensional gauge theory. 

As before, the D6-branes span (Euclidean) time and the spatial 123-directions in common 
with the D4-brane worldvolume, and lie at a fixed value of r. In parallel with the steps de- 
scribed by eqs. (2.11-2.14), we introduce isotropic coordinates {X,^2,f,<p} for the Euclidean 
black hole (6.2). Furthermore, we rescale to dimensionless coordinates 



deconf 



= MKK/27r . 



(6.5) 



A ^ C/tA , 



(6.6) 



in analogy with eq. (2.18). Now the embedding equation for r(A) becomes 




(6.7) 



14 



Note that the spectrum would again be calculated with Lorentzian signature. 
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Comparing this result with eq. (2.19), we see that the sign of the right-hand side has changed. 
As we will see below, this change of sign indicates that the brane now bends towards U = 11^: 
as expected, the black hole exerts an attractive force on the D6-brane, in contrast to the 
repulsive force in the D4-soliton background. 

If we look for asymptotically constant cmbcddings such that r(A) — > Too = L/Ut as 
A — >■ oo, we find the same long-distance behaviour as before: 

c 

r ~ roo + ■ (6.8) 

These asymptotic parameters are related to the quark mass and chiral condensate as before, 
except for the substitution of for U^^: 



(V'V) = -87r'^.ro6f/^c. (6.9) 



Again, the embedding equation can be solved numerically for any value of Too, with the 
asymptotic boundary conditions implicit in (6.8). Inspection of the results reveals that, just 
as in reference [14], two qualitatively different types of profile are possible: those for which 
the probe brane intersects the horizon, and those for which it does not. As we vary the 
boundary condition Too, we find that there is a unique embedding for sufficiently large or 
small values. In the case of small Too, the unique embedding is such that the D6-brane falls 
into the horizon, whereas in the large-roo case it does not, despite, of course, still being 
attracted by the black hole. For an intermediate range with roo ~ 1, we find that more than 
one embedding is possible. In particular, for a given Too, we find both cmbcddings which 
intersect the horizon and those where the D6-brane smoothly closes off before reaching the 
horizon. Furthermore, we also find that there are typically multiple solutions of at least 
one of these varieties of embedding; figure 6 displays some representative examples. These 
results imply that in varying from large to small Too, there will be a phase transition in the 
behaviour of the quarks. This transition will occur in the intermediate regime, where, for 
a fixed temperature and quark mass, the condensate in the gauge theory can have multiple 
values. Of course, the true ground-state condensate will be selected by minimizing the energy 
density. 

Combining various results — see eqs. (2.6), (6.3), (6.4) and (6.9) — the relation between 
the temperature and r^o is found to be 

T = ^, withM^.A!VV|K. (6.10) 

Here M can be interpreted, for a given quark mass, as the typical meson mass scale in 
the regime where supersymmetry is restored in the zero-temperature gauge theory — see the 
discussion below eq. (3.16). In the zero-temperature analysis of previous sections we regarded 
the variation of r^, as a variation of the quark mass, while [/kk or, equivalently, the KK scale, 
was fixed. In this section we wish to study the phase structure of the gauge theory as the 
temperature is varied, keeping the microscopic gauge theory parameters, and in particular 
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Figure 6: Some representative D6-brane embeddings from the region in which c is multi-valued. The 
three profiles in each figure have the same asymptotic value of ; in (a) two of them avoid the horizon 
whereas in (b) two fall into the horizon. 



the quark mass, fixed. For this reason, since in this phase the physics will only depend on 
the ratio T/M, as suggested by eq. (6.10), we will think of the variation of r^o as a variation 
of the temperature. 

The behaviour of the condensate (represented by c) as a function of the temperature is 
shown in figures 7a and 7b. These plots display some interesting features. First note that 
c — > both as T ^ and as T ^ oo. In the high-temperature regime, the unique embedding 
is such that the D6-brane falls into the horizon. By following the solid line in the figure, we see 
that these embeddings only exist down to some some minimum temperature, T ~ 0.9869M. 
Similarly there is a unique embedding at sufficiently low temperatures such that the D6-brane 
does not reach the horizon, and following the dashed line shows that embeddings of this form 
only exist up to T ~ 1.0202M. As shown most clearly in figure 7b, multiple solutions exist 
in the intermediate temperature range between the two values given above. 

Naturally, where more than one solution to the embedding equation exists, the ground- 
state profile will be that selected out by lowest energy density. Since we are considering static 
configurations, the energy density is given by 

S = - j dXd^zC = AttTj^^U^ j dA ^1 - [j^ j X^^/l + r'^ , (6.11) 

where A^iin is either zero, or the point where the brane meets the horizon, depending on 
which variety of embedding we are considering. This integral is divergent, so to compare the 
energy densities of different embeddings we regularize £ by using as a reference the embedding 
r(A) = 0, which runs from the horizon (at A = Ah, say) to infinity. Since this is a numerical 
calculation, we also introduce a cut-off, A^atch) in A at the point where, by imposing the 
appropriate boundary conditions, we match the numerical profile to the asymptotic solution. 
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Figure 7: The quark condensate (^'V') oc c (plots (a) and (b)), and the energy density £ (plots (c) 
and (d)) as functions of temperature. The solid (dashed) lines describe the family of embeddings that 
do (not) intersect the horizon. Plots (b) and (d) show in more detail the regions of plots (a) and (c) 
where the two branches intersect. The vertical line lies at T = Tfunj, the temperature at which the 
phase transition in the behaviour of the quarks takes place. 



r ~ Too + c/A . There is an infinite contribution to £ from the remainder of the profile, but 
when regularized the leading behaviour of this correction can be shown to be 



2A, 



+ c(A-Lj 



(6.12) 



match 
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Therefore, to be precise, we compare the energy densities by calculating 




(6.13) 



We then trust that this provides a valid basis for the comparison of different energy densities 
when the differences in £ are greater than 0(X^^^^). This is the case here, as we use A^atch = 
100. 

The energy density £ is plotted as a function of temperature in figures 7c and 7d. Starting 
from low temperatures, figure 7d shows that, as T is increased to that in the intermediate 
region, the physical embedding remains that which does not reach the horizon. At T ~ 
1.0202M, the energy density for this embedding and that which falls into the horizon match. 
For T > 1.0202M, the energetically favoured embedding is one which falls into the horizon. 
Therefore, at 

T = r,,„d ^ 1.0202M , (6.14) 

there is a phase transition in the behaviour of the fundamental quarks. This phase transition 
is signalled by a discontinuity in {tpip)(T), as shown in Figure 7b, as well as a discontinuity in 
the specific heat d£{T)/dT, as can be seen in Figure 7d. Note that, to within the precision 
of our numerical calculations, Tf^^^ cannot be distinguished from the maximum temperature 
for which a horizon-avoiding embedding is possible. 

From the string description with the D4 throat geometry, the basic physics behind this 
transition is straightforward. Increasing the temperature increases both the radial position 
of the horizon and energy density of the black hole. For small enough a temperature, the 
D6-brane is pulled towards the horizon but its tension balances the gravitational attraction. 
However, above the critical temperature Tfu^d, the gravitational force has increased to the 
point where it overcomes the tension and the branc is drawn into the horizon. 

Another feature that characterizes the present phase transition is that the meson spec- 
trum changes at T = Tf^^d- The discrete spectrum of mesons of the low-temperature phase 
becomes a continuum of excitations in the high-temperature phase where the D6-brane em- 
bedding extends to the horizon. This behaviour is similar to that of the glueball spectrum at 
the deconfinement phase transition in the gauge theory, discussed previously. Note, however, 
that for temperatures above Tdoconf we are, by assumption, in the deconfined phase of the 
gauge theory, and so it is possible to introduce free quarks; these correspond to open strings 
extending from the D6-brane down to the horizon. If Tjeconf < Tfund there is then a range 
of temperatures, T^econt < T < Tf^^^, in which there are free quarks but also a finite mass 
gap, mes{T) > 0, associated with the introduction of such quarks into the system. In this 
regime, there still exists a (discrete) spectrum of bound states with M < 2m^g. Note also 
that the deconfined adjoint degrees of freedom still generate a finite screening length beyond 
which the force between a quark and an anti-quark vanishes [26]; in the string description 
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this happens because the string that joins the quarks together snaps into two pieces, each 
of them going from the brane to the horizon, that can be separated with no cost in energy. 
Above Tfund, however, meff(T) = and the spectrum of fundamental matter field excitations 
becomes continuous. 

Up to this point, we have discussed two independent phase transitions: the deconfining 
phase transition at T = T^ccont, for which the holographic dual description involved a discrete 
change of the bulk spacetime geometry, and the quark phase transition at T = Tf^nd, for 
which the holographic description involved a discrete jump in the embedding geometry of 
the D6-brane. However, if the latter is to be realized as a separate phase transition, then it 
must be true that Tf^nd > ^deconf-^^ In terms of the microscopic parameters, this condition 
translates into 

"iq > ^5ymA^cMkk . (6.15) 

In the zero-temperature theory, this is roughly the quark mass scale at which supersymmetry 
is restored — see the discussion below eq. (3.16). 

If ^fund < Tdoconf ) thc rclevaut spacetime geometry at T = T(^^^ is the D4-soliton (6.1) 
and the ambiguity in the D6-brane embedding discussed above is not relevant. Alternatively, 
the gauge theory is still in a confined phase and so the chemical potential to introduce a 
free quark remains infinite; in the dual description this is because, since there is no horizon, 
an open string must have both ends on the D6-brane, and so represents a quark/anti-quark 
bound state. In this case, the two phase transitions are realized simultaneously at T = T^cconf- 
An important example of this behaviour is when = 0. In this case, the low-temperature 
phase is characterized by spontaneous chiral symmetry breaking, as in section 2. In the 
high-temperature phase, mq — is equivalent to T — oo (see eq. (6.10)), in which limit 
the chiral condensate vanishes (see figure 7a). Hence the deconfining phase transition is also 
accompanied by the restoration of chiral symmetry in this model. 

To close the section, we wish to stress to the reader that while these phase transitions 
are of inherent interest in understanding the model under study here, they take place in 
a temperature regime where the gauge theory is intrinsically five-dimensional, rather than 
behaving in a four-dimensional way, as can be seen by comparing the transition temperatures 
given in eqs. (6.5) and (6.15) to the compactification scale Mkk- 

7. D4/D6/D6-intersections 

In the discussion of embeddings in section 2, we noted that the solution r(A) = of eq. (2.19) 
is not a physical one (for J7kk / 0) if it terminates at the origin of the C/r-plane because 
it would correspond to an open D6-brane. However, this solution can be extended through 
this point to construct a physically acceptable embedding. To see this, recall that the full 

^^Note that, in this regime, the fundamental fields still 'feel' the deconfining transition at T = Tdeconf- For 
example, there is a discontinuity in the chiral condensate {'ijjtp), which jumps from a negative to a positive 
value at this phase transition. 
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D6-brane embedding is further specified by the conditions <p = (po and r = tq, where 4>o and 
To are arbitrary constants. We can therefore join together a solution with r = tq to another 
one with t = tq + St/2 (see eq. (2.3)) to obtain a complete regular solution that describes a 
D6-brane whose projection on the Ur-plane is a straight line that passes through the origin 
and intersects the boundary at two antipodal values of r. This means that this solution 
describes an intersection of D4-branes with two D6-branes, or, more precisely, with one 
D6-brane and one anti-De-brane.-*^^ Note that these embeddings bear a striking resemblance 
to those discussed in [13]. The parallel includes the facts that both constructions involve two 
boundary defects and both preserve the chiral symmetry. 

To better understand the defect theory corresponding to these D4/D6/D6-intersections, 
we consider a family of more general embeddings that describe D6-branes that lie at r = 
(or alternatively, = = z^) and stretch along some curve U = U{t) in the J7r-plane, 
intersecting the boundary at two, not necessarily antipodal, values of r. Such D6-branes wrap 
a maximal S'^ within the in the metric (2.1), so the induced metric on their worldvolume 
is 
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Since the Lagrangian density on the D6-brane 
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(where U' = dU/dr) does not depend explicitly on r, it follows that 
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is also r-independent. Let Uq > C/kk be the minimum value of U{t). By symmetry we can 
assume that this value is reached at r = 0, that is, that we have U{0) = Uq and U'{Q) = 0. 
Integrating the equation above we then find 



f{u) = lhy^ g'/\h) £dx 
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where we have introduced 
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(7.5) 



Note that the rescaled angle f has period 2-k. Eq. (7.4) specifies the D6-brane profile in 
the J7r-plane. The D6-brane intersects the boundary at f±{b) = ibf(n —>■ 00). These two 

^^Technically, the reason is that the orientations on the intersections of the D6-brane with the boundary at 
T = To and T = To + 5t/2, induced from a given orientation on the D6-brane, are opposite to each other. 
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points are implicitly fixed by Uq, the 'lowest' point of the D6-brane in the t/r-plane. In 
particular, their separation is a monotonically decreasing function of b. It is easy to show 
that f±{b — 1) = ±7r/2, so the two points become antipodal as the lowest point of the D6- 
brane approaches the origin of the U r-plane. In this limit the projection of the D6-brane on 
this plane is just a straight line, as discussed above. In the opposite limit in which the (lowest 
point of) D6-brane approaches the boundary, the separation between the intersection points 
becomes smaller and smaller. In the strict limit b ^ oo the two coincide and the D6-brane 
disappears. 

Now, these embeddings seem to portend fascinating physics for the dual five-dimensional 
gauge theory coupled to two defects or, more precisely, to a defect and an anti-defect. It 
appears, however, that this physics is intrinsically five-dimensional. Since the focus of the 
present paper is the four-dimensional QCD-like physics of the original model, we have not 
fully explored this system. However, we offer some preliminary remarks in the following. 

The deflection of the D6-brane in the r-direction, which corresponds to a worldvolume 
direction of the D4-branes, is reminiscent of certain D5-brane embeddings in AdS^ studied 
in ref. [27]. These authors found (non-supersymmetric) embeddings in which the D5-brane 
began at large radius running along r, but were deflected to run parallel to the AdS^ horizon at 
some finite radius. At large r, these configurations can be understood as D3/D5 intersections 
that are T-dual to that represented in (1.1). The dual gauge theory involved coupling Af = 4 
SYM to a (2 -I- l)-dimensional defect which broke all super symmetries. The deflection was 
interpreted as giving an expectation value to an operator in the defect theory. 

In the present case, the D4-soliton background breaks the supcrsymmetry but with both 
defects, supcrsymmetry would still not be restored in the limit Ukk 0;^'' one way to under- 
stand this in the field theory is to note that the half of the supersymmetries of the ambient 
theory preserved (broken) by the defect are precisely those that are broken (preserved) by 
the anti-defect.^^ The embeddings described by eq. (7.4) would correspond to a state where a 
symmetric pair of operators in each defect theory simultaneously acquire identical expectation 
values. Indirectly, the D6-branc embedding would also induce different expectation values of 
gauge theory operators, e.g., Tr F^, in the two regions on the r-circlc divided by the defects. 

At the microscopic level, there are two sets of fundamental matter fields in the field 
theory, one associated with each defect. This is realised in the string description by the fact 
that, despite there being a single D6-brane, its intersection with a constant-t/" > Uq slice 
consists of two disconnected pieces, corresponding to two independent sets of active degrees 
of freedom at the corresponding energy scale. These degrees of freedom seem to disappear 
for energy scales associated to radial positions U < J7o) as in [27]. As discussed above, the 
D6-brane disappears in the limit f_|_ — f_ — > 0. Consistently, in this limit Uq ^ oo and the 
degrees of freedom above are completely removed. In the gauge theory, this corresponds to 
the limit in which the defect and the anti-defect are brought on top of each other. 

^^Note that the supersymmetry breaking of the background is not essential to constructing the general 

embeddings, i.e., one can simply replace g{x) by 1 everywhere in eq. (7.4). 

^^Defect/anti-defect configurations, and more general multiple defect configurations, were discussed in [28]. 
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This situation may be compared with that where both defects are associated with a 
'standard' D6-brane embedding as described in section 2; one at f = f+ and the other at 
f = f_. Of course, in this case there are also two sets of independent degrees of freedom 
associated with open string excitations on both the D6- and D6-brane.^^ However, one 
difference is that the spectrum of these excitations will now be independent of the defect/anti- 
defect separation, f+ — f-?^ In addition, one would also find excitations corresponding to 
open strings stretching between the two separate branes. In contrast, in the case of a single 
D6-brane, the two defects can interact through open string modes that propagate along the 
brane from one defect to the other. 

Of course, one should ask the question which embedding minimizes the energy density 
of the system. That is, does energetics favour the smooth embedding which joins the two 
boundary defects or that in which there are two independent branes? As a preliminary step 
in this direction, we compared the energy density for the 'standard' embeddings of the D6- 
and D6-brane pair with Too = to that of the smooth embedding above with Uq = [/kk (and 
hence f + — f_ = tt) . Our numerical results indicate that the smooth embedding produces a 
lower energy density and hence describes the true vacuum configuration. Note then that with 
both defects corresponding to the D6- and D6-branes, the vacuum preserves chiral symmetry. 

Presumably if Uq is increased above Ukk, the smooth embeddings described above would 
continue to be energetically favored as it seems the energy density for these solutions must 
decrease as f+ — f_ decreases while that for the standard embeddings remains unchanged. 
Note, however, that the above are just two particular cases of a more general family of 
embeddings in which the brane bends simultaneously in the 89- and the ?7r-planes, that is, 
embeddings with nontrivial functions r = r(A) and r = t(A) (using the notation of section 
2). In general, we expect the minimum-energy two-defect embedding is of this type. 

These defect theories might be explored in other ways as well. For example, one might 
also consider varying the parameter for the two defects, which may be done independently 
for each of the defects. Finally, of course, one might also consider the effects of a finite 
temperature, as in section 6. 

8. Discussion 

We have analyzed in detail the case of a single D6-brane, which corresponds to a single 
flavour in the gauge theory. In this case the gauge theory enjoys (at large N^) an exact 
U{l)y X ?7(1)a symmetry. The first factor is just the flavour symmetry, which in the string 
description corresponds to the gauge symmetry on the worldvolume of the D6-brane. The 
second factor is a chiral symmetry that rotates the left- and right-handed fermions, and 
V'r, with opposite phases; it also rotates by a phase the adjoint scalar X = + iX^, as is 

^^Note that in this case one may also reaUse a defect/defect (as opposed to a defect/anti-defect) theory by 
considering two D6-branes. 

^"in the approximation in which interactions between the branes, which are suppressed in the large- iVc limit, 
are neglected. 
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clear from the fact that in the string description the U{1)a symmetry is just the rotational 
symmetry in the 89-plane. 

Like in QCD, we expect the U{1)a symmetry to be spontaneously broken by a chiral 
condensate, (ipip) ^ 0, and hence the existence of the corresponding Goldstone boson in the 
spectrum; this is the analog of the QCD rf particle (which is masslcss in the ^ oo limit), 
but, in an abuse of language, we have referred to it as a 'pion'. The string description indeed 
confirms the field theory expectation: for zero quark mass, there is precisely one massless 
pseudo-scalar (in the four-dimensional sense) fluctuation of the D6-brane. If mq > this pion 
becomes a pseudo-Goldstone boson. Its mass, as well as the chiral condensate and the pion 
decay constant, can be computed in the string description, and we have shown that these 
three quantities and the quark mass obey the Gell-Mann-Oakcs-Rcnner relation (1.2) [18]. 
Although this is a reassuring result, one should keep in mind that, just like in field theory, it 
is essentially a consequence of the symmetry-hreakmg pattern. 

In the case Nf > 1 we provided a holographic version of the Vafa-Witten theorem, 
which states that, in a vector-like theory with zero 0-angle (like QCD) and Nf quark flavours 
of identical masses > 0, the global, U{N{) flavour symmetry cannot be spontaneously 
broken. As we discussed, it is a priori unclear, on field theory grounds, that the Vafa-Witten 
theorem does apply to the actual gauge theory on the D4-branes, because of the presence 
of Yukawa-like couplings to scalar fields of masses Mkk ~ Aqcd • The result from the string 
analysis therefore established that the presence of these fields does not alter the vacuum 
structure of the theory to the extent of violating the Vafa-Witten theorem. 

The global U{Nf) symmetry of the gauge theory becomes the gauge symmetry on the 
worldvolume of Nf D6-branes in the string description. The holographic version of the Vafa- 
Witten theorem is the statement that the minimum-energy configuration for the A^f D6-branes 
is a [/ ( Af)-preserving one in which all branes lie on top of each other. It should be emphasised 
that, unlike the Gell-Mann-Oakes-Renner relation, the Vafa-Witten theorem is not based on 
a symmetry (-breaking) argument, but depends on the dynamics of the theory. In the field 
theory, the proof relies on the reality and positivity of the fermion determinant that arises in 
the path integral formulation; in the holographic version, it relies on details of the dynamics 
of multiple D6-branes. 

One result of our analysis of the Nf > 1 case is the existence of N^ independent, massless 
fluctuations of the D6-branes if = 0. It is tempting to interpret the associated N^ 
massless pseudo-scalar particles in the spectrum of the gauge theory as the Goldstone bosons 
of a spontaneously broken U{Nf)f^ symmetry, to which the t/(l)A symmetry would putatively 
be enhanced for Af > 1. However, only one of the A^^ states above is a true Goldstone boson, 
for two reasons. First, the five-dimensional theory contains a term of the form ^iXij)^, where 
i = l,...,Af is the flavour index. Since X only transforms under the U{1)a subgroup of 
U{Nf)A, the theory is only invariant under U{1)a x U{Nf)Y, which is spontaneously broken to 
U{N()y. Although A is a massive field that can be integrated out, its mass Mx ~ ^kk ~ Aqcd 
is of the same order as the strong coupling scale, so its integration should not lead to a real 
suppression of the symmetry-breaking operator above. 
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The second reason that makes implausible the interpretation of all but one of the massless 
modes above as Goldstone bosons is that they have non-derivative interactions. Indeed, the 
effective, interacting Lagrangian for these modes is obtained by expanding the non-Abelian 
Dirac— Born-Infeld action for the Nf D6-branes [23] beyond qiiadratic order and then dimen- 
sionally reducing it to four dimensions, as we did in Section 4 in the Ni = 1 case. Since 
the action for the D6-branes contains non-derivative interactions, such as commutator terms 
squared, so does the effective four-dimensional Lagrangian. These terms are only absent for 
the mode that multiplies the identity generator of U {Nf) . The caveat as to why this argument 
makes the Goldstone boson interpretation for the other modes implausible, but not strictly 
impossible, is that the coefficients in the effective Lagrangian might conspire so that, once all 
contributing diagrams are included, only momentum-dependent parts survive in any scatter- 
ing amplitude (as in the linear sigma-model for the pion Lagrangian [29]). This possibility 
seems rather unlikely. 

Thus we conclude that the string description predicts the existence of massless pseudo- 
scalar particles in the spectrum of the quantum gauge theory in the large- AT^ limit. Once sub- 
leading effects in the 1 /N^ expansion are included, we expect these particles to acquire a small 
mass. The U{l)j^ symmetry is anomalous at finite N^, so its associated Goldstone boson is not 
exactly massless at finite N^; presumably, this anomaly can be analysed in the supergravity 
description along the lines of [6, 30, 31, 32]. The other — 1 modes are likely to acquire a 
mass due to closed string interactions between the D6-branes, which are a ~ 1/-/Vc effect. 
Prom the gauge theory viewpoint, we have found no obvious reason for the masslessness (or 
lightness, at finite N^) of these — 1 pseudo-scalars, since we have argued that they are not 
Goldstone bosons. ^""^ Recall that, from the string viewpoint, the reason for this is that the only 
difference between the non-Abelian and the Abelian DBI actions consists exclusively of terms 
that involve commutators of fields. This can be regarded as a consequence of the fact that 
the seven-dimensional D6-brane action is T-dual to the ten-dimensional action of D9-branes; 
all commutators of scalar fields in the D6-branc action originate from commutators of gauge 
fields in the non-Abelian field strength on the D9-branes. Thus, it is ien-dimensional U{N{) 
gauge invariancc that seems to be responsible, from the string viewpoint, for the lightness of 
the gauge theory scalars.^^ Given the importance in high energy physics of mechanisms that 
naturally allow small scalar masses, it would be interesting to understand how to rephrase the 
constraints imposed by ten-dimensional U{N{) gauge invariance purely in four-dimensional 
field theory terms. 
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